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Abstract 

We calculate the shear (rj) and bulk [Q viscosities of a weakly coupled quark gluon plasma at 

o , 

CO the leading-log order with finite temperature T and quark chemical potential ji. We find that the 

in ' 

■ shear viscosity to entropy density ratio n/s increases monotonically with ii and eventually scales 

■ as (/z/T) 2 at large \i. In contrary, Q/s is insensitive to \i. Both rj/s and Q/s are monotonically 
decreasing functions of the quark flavor number Nf when Nf > 2. This property is also observed 
in pion gas systems. Our perturbative calculation suggests that QCD becomes the most perfect 

H ■ 

(i.e. with the smallest rj/s) at ji = and Nf = 16 (the maximum Nf with asymptotic freedom). 
It would be interesting to test whether the currently smallest n/s computed close to the phase 
transition with /x = and Nf = can be further reduced by increasing Nf. 
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I. INTRODUCTION 



Viscosity, diffusivity, and conductivity are transport coefficients which characterize the 
dynamics of long wavelength and low frequency fluctuations in a medium. The quantity shear 
viscosity (77) per entropy density (s) has attracted a lot of attention because of the intriguing 
conjecture that rj/s has a minimum bound of 1/4tt for all systems This conjecture is 
inspired by the anti-de-Sitter space/conformal field theory correspondence (AdS/CFT) (2 - 4 1 
which is rooted in string theory. Surp risingly, the hot and dense matter produced at RHIC 
5-8] (for reviews, see e.g. Ref. [9 12]) just above the phase transition temperature (T c ) has 



rj/s = 0.1 ± 0.1(theory)±0.08(experiment) [13j, a value close to the conjectured bound. A 
robust limit of 1/(4-71") < rj/s < 2.5/(47r) at T c < T < 2T C was recently extracted from a 
VISHNU hybrid model jl^] and a lattice computation of gluon plasma yields rj/s — 0.102(56) 
at temperature T = 1.24T C jl5 ]. 

The QCD transport coefficients have also been studied in other temperatures. When 
T ^> T c , the t] of a weakly interacting quark gluon plasma is inversely proportional to the 
scattering rate, r\ oc 1/T oc l/a^lna" 1 [161] . where a s is the strong coupling constant. The 
bulk viscosity ( is suppressed by an additional factor of (T^) 2 , arising from the response of 
the trace of the energy momentum tensor (T^J to a uniform expansion. Thus, ( vanishes 
when the system is "conformal" or scale invariant. For a gluon plasma, the running of the 
coupling constant breaks the scale invariance. Thus, oc (3 (a s ) oc a 2 , ( oc a 2 / In a~ l 
171 ] . In the perturbative region, (/r] oc a; 4 <C 1. When T <C T c , the effective degrees of 
freedom are pions. In the chiral limit (u and d quarks are massless), rj/s oc / 4 /T 4 18] and 
(/s oc T 4 // 4 [19| where is the pion decay constant. A compilation of perturbative QCD 



calculation of rj and ( can be found, e.g., in Ref. |20|, |21| . Most of these calculations are 
performed with finite T but zero quark chemical potential /i. 

The purpose of this work is to extend the previous perturbative QCD calculation of rj 
and ( to finite \i at the leading-log (LL) order. At this order, we find 77 oc l/a^lna^ 1 
and ( oc c^/lna^ 1 , same as in the limit of /1 = 0, which give parametrically the dominant 
contribution at asymptotically large (compared with the Aqcd scale where QCD becomes 
non-perturbative) T or /i. The vacuum in our calculation has no spontaneous symmetry 
breaking, thus, it cannot be applied to the color superconducting phase in the fi/T — > 00 
limit. In the context of finding the minimal rj/s and hence the most "perfect" fluid, we 
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explore whether rj/s can be further reduced by varying \x and the quark flavor number Nf 
in the hope that our perturbative calculation can shed light on the non-perturbative region 
near T c where QCD is found to be the most perfect matter ever produced in Nature. 



II. EFFECTIVE KINETIC THEORY 



r) and ( can be calculated through the linearized response function of thermal equilibrium 
states using the Kubo formula. In the leading order (LO) expansion in the coupling constant, 
the computation involves an infinite number of diagrams 22M25I] . However, it have been 
shown thatthe —ion of the LO digrams in a weakly eonpled <t theory HQ or in 
hot QED |27j] is equivalent to solving the linearized Boltzmann equation with temperature- 
dependent particle masses and scattering amplitudes. This conclusion is expected to hold 
in perturbative QCD as well. 

The Boltzmann equation of a quark gluon plasma describes the evolution of the color 
and spin averaged distribution function f£(x) for particle a: 



dt 



(1) 



where f£(x) is a function of space-time x M = (t, x) and momentum p M = (E p , p). For the LL 
calculation, we only need to consider two particle scattering processes denoted as ah -H- cd 
with the collision terms in the form 



Cabled — / dTab -^. c d 

J k]_k2k3 

where F 9 = 1 + f 9 and F 9 ® = 1 - /«® and 

3 
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(2) 



dT 



1 



ab—tcd 



2E n 



\M, 



ab^cd\ 



n 



d 3 L 



(3) 



A (2tt) 3 2£, 

where \M a ^ C( i\ 2 is the matrix element squared with all colors and helicities of the initial 
and final states summed over. They are tabulated in Table [J in the Appendix. Then the 
collision term for a quark of flavor a is 



NX, 



1 



'A''/'. " c^Qaqa^qaqa ^ q a q a -i^,q a q a + ^' ggJ^q a q a + C ^q a g-^-q a g 

^ ] (Cq a q b 4^,q a q b + C ' q a q b ^q a q b + ^ q b q b -^q a qa ) ) 
b,by^a 



(4) 



where N q = 2 x 3 = 6 is the quark helicity and color degeneracy factor and the factor 1/2 
is included when the initial state is formed by two identical particles. Similarly, 



NgCg — yCgg^gg + ^J(C ' gq a ^gq a + Cgq a ^gq a + Cg a g a 



(5) 



where N g = 2 x 8 = 16 is the gluon helicity and color degeneracy factor. In equilibrium, the 
distributions are denoted as f q ^ and f 9 , with 

1 



fa 



e u-p/T _ l 
1 



(6) 
(7) 



e (u-Ffn)/T _|_ X ' 

where T is the temperature, u is the fluid four velocity and \i is the quark chemical potential. 
They are all space time dependent. 

The thermal masses of gluon and quark or anti-quark for external states (the asymptotic 



masses) can be computed via [28 



m„ 



2g 2 r d?p 
~cU J (2vr) 3 2E p 



NgCjj; + N f N 9 C F (J* + f} 



m 2 g = 2C F g 2 



d 3 p 



(2tt) 3 2E p 

where = 8, Ca = 3, and Cp = 4/3. This yields 



(2/* + r v + /«; 



(8) 
(9) 



C 



m 



A „2rp2 



N f Cf 



rn 



6 y 16 y v 



7T^ 



;C F g 2 T 2 + 



7T- 



(10) 

(11) 



where we have set E p = |p| in the integrals on the right hand sides of Eqs. dHJ and (Q. The 
difference from non- vanishing masses is of higher order. 



A. Linearized Boltzmann Equation 

To extract transport coefficients, it is sufficient to consider infinitesimal perturbations 
away from equilibrium which have infinite wave lengths. Using the Chapman- Enskog ex- 
pansion we linearize the Boltzmann equation to the first order in the derivative expansion in 
x. Thus, we only need the thermal equilibrium distributions f 9 and f g,q on the left-hand-side 
of the Boltzmann equation. We can also make use of the zeroth-order energy-momentum 
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conservation relation, d tl T^ tlv = 0, to replace the time derivatives dT/dt and dfi/dt with 
spatial gradients: 



dT m fdP\ _ 



7ft 



&p\ fdp\ 

» I "57 + 77T 



\ dn ) 



Vu, 



(12) 



where n 



n 5 is the baryon number density. We work in the local rest frame of 



the fluid element where u = (1,0,0,0) which implies d^u = after taking a derivative 



on ^(xju^x 



1. For the right-hand- side of the Boltzmann equation, we expand the 



distribution function of particle a as a local equilibrium distribution plus a correction 



and x a can De parametrized as 



A a ( p ) B a (p)^^^ 
Villi H ^—P[iPj\V\iUj] 



(13) 



(14) 



rj-\ <> <> 1 rj-\ 

where i,j = 1,2,3 and [...] means the enclosed indices are made symmetric and traceless. 

Applying Eq.ffTBl to the right-hand-side of the Boltzmann equation and equating it to 
left-hand-side, we get linear equations for B a (p) and A a (p). For f 9 , we obtain 

N f 



P[iPj) 



E g 1 

fBFBNg 



1 

2 99*^99 ^ Z-^i \ ma*+gq a ~ 9q a ^91 a ^ ^QaQa^gg J J > \ LU ) 
a=l 



where 



B ab^cd(P) 



dY ab ^ cd r k J b k2 F c p F* 3 [-B^h) - Sj(fe) + B§(p) + Sj(fcs)] • (16) 



with Bf-ip) = B a (p)p\iPj] and we have suppressed the p dependence in 5?. and 5 a when 
there is no ambiguity. Similarly, for f q , we obtain 

1 



P[iPj] 



E « 1 



+ r'j_ _ — i— :r«j _ + B tJ 
fqpq m |_2 qiq^qiqi qiqi^qiqi ~ o 99**q\q\ ~ qig^qig 

J Q — 



N f 



~^ i^qiqa-n-qiqa ^qiqa-^qiqa + B q a qa<^qiqi) 



a=2 



(17) 



The corresponding equation for can be obtained from the above equation by interchanging 
q -H- g. The linear equations for £? a (p) is relevant to the shear viscosity computation. They 
can be written in a compact form 

\S ri )=C ri \B), (18) 



with P[iPj] taken as the source IS"?,) for the shear viscosity. 

Similarly, the linear equations for A a (p) for quark and gluon are given by 



V 
~3 



,<9m 2 



I T7i rp2 ^ ,,u 9 2 9 



<9m 2 \ /&p\ 



— a, 



<9P 

cm 



1 \ 



0=1 



(19) 



and 



V 
~3 



_ , E -2 _ T2 dm 2 q 2 dm 2 q \ ( dP\ 



Q T 2 



W) \de) n + E * 



d_Eg 
djj, 



— a, 



dP\ 
dn J 



Eg 1 



fiFi N q 

~l~ ^ ] (^-giga^giga "I" ^<?i<ja-H>gi<?a "I" Ag a q a 



a=2 



(20) 



where a 5 = 0, a 9 /g = ±1 and 

A^M = I dT ab ^ cd fZJ b k2 FlF p d [-A a (h) - A b (k 2 ) + A c (k 3 ) + A d (p)] . (21) 

These equations are relevant for the bulk viscosity computation. They can be written 
compactly as 

\S C )=C ( \A). (22) 

Using the expression of the thermal mass in Eq. (jHJ), we can evaluate the prefactor of 
(dP/de) n in Eq. (JE) as 

3 



L- rrV — 'g n 2 9m l 



P 2 - P(g' 



2 i ~ 2 

P +m , 



C A T 2 , N f t F 

6 6 



T 2 + — fj 2 



7T Z 



where a s = g 2 /4:ir and the QCD beta function for g(jl= a/T 2 + y^ 2 ) is 



- 4 / 2iV / -33 



r,2 



16tt 2 



(23) 



(24) 



Similarly we obtain the prefactor of (dP/de) n in Eq. (120 j) as 



dT 2 



dfi 2 



P 2 - /% 2 ) 



2 , ~ 2 

P +«V 



(25) 



As will be shown later, (dP/de) n — 1/3 oc /3(g 2 ) and (dP/dn) e oc (3(g 2 ), thus the source 
of the bulk viscosity oc (3(g 2 ). This means |S^) is, as expected, proportional to the 
conformal symmetry breaking. 
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B. Energy-Momentum Tensor and Quark Number Density 



difference in metric with 

rpflU 



The energy-momentum tensor of the kinetic theory can be written as [22| (note the sign 

3) 

d 3 p f a (p,x) 



(2tt) 3 E a 



(26) 



where a sums over the gluons and Nf flavors of quarks and anti-quarks. This equation 
expresses the total energy momentum tensor of the system as the sum of individual quasi- 
particles. There is no u^u v term on the right hand side because energy momentum conser- 
vation cannot be satisfied unless this term vanishes. In principle, one expects the form of 
Eq. (|26|) (and kinetic theory itself) will be no longer valid at higher orders in the expansion 
of coupling constant, however. Eq. (126 j) does reproduce all the thermal dynamical quantities 
of QCD at 0(a s ) correctly [3o|. 

Expanding T^ v to the first order of x\ we have the fairly simple expression 

/r/ 3 w / BE 2 8E 2 \ 

__f_ rFV ^ _ uVT 2 ^ - uVy^f) . (27) 

In deriving the above equation one needs to carefully keep track of the implicit distribution 
function dependence in E a through the thermal mass This expression can then be 
matched to hydrodynamics. 

In hydrodynamics, T^ v depends on T, /i, and u^. In the local rest frame of the fluid 
element with u(x) = 0, the most general form of T u at the order of one space-time derivative 
is (assuming parity and time reversal symmetry) can be decomposed into shear and bulk 
viscosity terms 

ST« = -277 2 ~ • u) - tf y V ■ u. (28) 

(Recall that dT/dt and dfi/dt can be replaced by V • u as shown in Eq. (I12j) .) Then using 
Eq. ffMl) for x an d comparing Eqs. fl27l) and fT28|) . we obtain 

''-m^ N 'Iwk rF °^' Vih (29) 

and 

<^Ejj$wS f ° F ' A '- (30) 
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They can be written schematically as 

V=(B\S V ), C = (A\S' C ). (31) 

Note that \S' ( ) ^ \S C ), but we will show ( = (A\S' ( ) = (A|£ c ) later. 
We also need the total quark number density 

Expanding n to the first order of x"? we have 

^X>/|^"(f -«•)■ (33) 
Eqs. ( 1271 f33l) will be used in computation of bulk viscosity. 

C. Shear Viscosity 

It is well known that r\ should not be negative such that the second law of thermodynamics 
(entropy cannot decrease in time) is satisfied. This requirement is fulfilled by rewriting Eqs. 
PEED as 

V = (B\C V \B) (34) 

and showing that rj is quadratic in \B) with a positive prefactor such that it is bounded 
from below by zero. Indeed, these conditions are satisfied in our expression: 

1 aa^gg ' 2 ^ <2<>?6->3a?& <Za96-*<?a?& ~ r oiy gaqb-*qaqb "r" 0±y qag a ^gbqb) a ^b 

TJ, J y-^qq-^qq ' ^ qq-^rqq ~ qq-^rqq T gg^qq T ^±JS qg ^ qg T qg^qg ) i \°°) 

where there is no summation over a and b and we can just take (a, b) = (1, 2) and where 

D v aM = J 0^dr ab ^ cd fZJ b k2 FlF p d [B^h) + 2% (fa) - B^ih) - Bf {p)] 2 . (36) 

Once 77 has the standard quadratic form in \B), and 77 = (B\C V \B) = (B\S V ), we can use the 
standard algorithm to systematically approach 77 from below 20] . 



S 



D. Bulk Viscosity and the Landau-Lifshitz Condition 



For bulk viscosity, the collision kernel Q in Eq. (|22|) has two zero modes Ae and A n 
which satisfy 

C ( \A E ) =C c \A n ) = 0. (37) 
Ae arises from energy conservation 

A a E (p) = E a , a = g,q,q, (38) 

while A n arises from quark number conservation 

i4»(p) = 0, A*(p) = l, AUp) = -1. (39) 
We can use the Landau-Lifshitz condition 

5T 00 = (40) 

and 

5n = (41) 

to rewrite Eq. (130]) by adding linear combinations of 5T 00 and Sn: 

Na f d*p , a „ aAa \p 2 fdP\ ^OEl ,8El 



+ (?*L\ Ea 



~ I -57) \K-T^-^ 



3 \de J n \ a dT 2 ^ dfi 2 
= (A\S C ). (42) 
Substituting \S^) with Eq. ( )2"2"|) . ( becomes quadratic in \A) with a positive prefactor: 
C = (A\C C \A) 



99^99 ' 2 ^ fcfl^Wt ^ 3a3i,->9a3f> 0± ^q a qb^qaqb 0±y qaqa->qbQb j , h 



Nf(N f - 1) / c , 

+iVf ^-D^gq + D$q-+qq + 4:D^^ qq + ^D^gg-^qq + ^qg^qg + ^qg-^-qg J > (43) 



where 

1 /" <i 3 p 



= ^ / 7^d^c d fU b k2 F^F p d [A"(h) + A\k 2 ) - A%k 3 ) - A d (p)] 2 . (44) 
(dP/de) n and (dP/dn) e in Eq. (|42|) can be obtained from thermodynamic relations 

/dP\ nTP^ T + {nii-sT)P^ T -siiP^ 

U J . = ^ ' (45) 



9 



where the pressure P (T, //) can be read of from Eq. (1261) and where P,x,y = d 2 P/dXdY , 
n = aP/«9/i, s = <9P/<9T, e = -P + Ts + fin and CV = T(P r>T - P 2 iT /P MlM ). We have also 
used 



dP 
Ik 



(dP dP \ 

w 5T+ dji 5f * 

\ a j ofi / (5ri= |n (5T+ |2i jAt=0 



(46) 



Then it can be shown that 

(A E \S C ) = 0, (A n |<S c ) = 0. (47) 

Thus one can solve ( = (A\CAA) = (A\S^) for ( using the standard algorithm to systemati- 
cally approach ( from below 20 1. 



E. Numerical Results 



As mentioned in the introduction, the LL result has rj oc g~ A In -1 (1/g) and £ (X 
g 4 ln _1 (l/g). When fi/T <C 1, rj, £, and s all scale as T 3 from dimensional analysis. For 
T//i < 1, 77, C, and s should be even functions of fi because they should be invariant under 
fi — > —fi, i.e. the exchange of quarks and antiquarks. 

We first show the result for the entropy density s. Only the leading order s for free 
particles is needed: 

s = N gSg + NfN q ( Sq + s q ) , (48) 



where 



v l 



f3E p 



In (l 



-0E P 



- + In (1 + e 



(49) 
(50) 



g/3(E p TM) _|_ 

In Fig. [U s/T 3 is shown as a function of (fi/T) 2 . When ju/T< 1, s scales as T 3 and when 
T//i < 1, s scales as /i 2 T. This agrees with the expectation that s = when T = 0, and s 
increases with T for fixed fi. The entropy density s also increases monotonically with the 
number of flavors Nf. We stop at Nf = 16, just before the asymptotic freedom of QCD is 
lost when N f > 33/2. 

In Fig. [2j we show the normalized shear viscosity fj = (t] /T 3 )g 4 In g' 1 as functions of 
(fi/T) 2 for different Nf. When fi/T — > 0, 77 scales as T 3 from dimensional analysis. When 
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FIG. 1: (color online) s/T 3 as functions of (li/T) 2 for different Nt 
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FIG. 2: (color online) (r]/T 3 )g 4 In 5 1 as functions of (/i/T) 2 for different A^y 



T/ [a —} 0, rj cannot scale as /1 3 because it is an even function of li. Instead, 77 scales as \x /T. 
Technically, this is because f q F q oc 5 ((E — li) /T) = T5 (E — li) as T/li — > while the anti- 
quark and gluon contributions vanish. Thus, Eqs. fl29l) and fl35l) are solved with oc 1/T 
and 77 oc 1/T. Physically, this 1/T behavior emerges because 77 is inversely proportional to 
the collision rate which vanishes at T = 0. Also, 77 is monotonically increasing with Nf 
because the averaged coupling between gluons is stronger than those with quarks involved. 
Thus, the effective collision rate is smaller with higher li and higher Nf. 
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Around \i = 0, we make a Taylor expansion fj = a v + 6^(/x/T) 2 + • • • , where a v is f/ at 
zero quark chemical potential. The values of a v and b v for various Nf is tabulated in Table 
HTl Our a v is identical to AMY's to at least the second decimal place for all Nf computed in 
Ref. [l6]. The agreement is better than 1%. 

In Fig. [31 (f] / s) g 4 In g -1 is shown as functions of (fi/T) 2 . For a given coupling, the LL 
value of rj/s is the smallest at /x = 0, i.e. the fluid is the most perfect at \x = 0. This 



perturbative QCD result at high T is consistent with the observation of Ref. [3l| in the 
hadronic phase at low T. Thus, we speculate that this property might also be true near the 
phase transition temperature T c such that QCD has its local minimum, perhaps its absolute 
minimum as well, at T c with zero quark chemical potential. 

If the coupling g is held fixed, our 7]/s is monotonically decreasing with Nf for Nf > 2 but 
not for Nf = and 1 (there is a crossing between the rj/s of Nf = 1 and 2 at /i 2 /T 2 ~ 1.8). 
This pattern looks random, but interestingly, it is qualitatively consistent with the pion gas 
result of Ref. [18[ which has r]/s oc f 4 /N 2 T A oc N 2 /N 2 (we have used f 2 oc N c ). Again, 
this suggests that there is a natural connection between rj/s above and below the phase 
transition. This pion gas analogy can also explain why Nf = and 1 are special — there is 
no pion in these two cases. 

We also observe that when Nf > 8, the fluid can be more perfect than that of Nf = 0. 
It would be interesting to test this in lattice QCD to find a more perfect fluid than the 
currently evaluated Nf = case. It would also be interesting to investigate how the above 
qualitative Nf scaling changes due to the possible infrared fixed point for Nf > 12 where 
chiral symmetry is not supposed to be broken and hence no pions exist anymore (see, e.g. 



32J and references therein). 

In Fig. 1 we show the normalized bulk viscosity £ = (£/a 2 T 3 ) hig -1 as functions of 
(/i/T) 2 for different Nf. When fi/T — > 0, £ scales as T 3 from dimensional analysis. When 
T I \i — > 0, ( scales as /i 2 T. Technically, this is because as T — > 0, f 9 = p = 0, the 
dominant contribution in Eq.( j43l) comes from the scattering between quarks. In quark 
scattering, the combination K = fkJk 2 Fk 3 F p [A a (h) + A\k 2 ) - A c (k 3 ) - A d (p)] 2 = 0(T 2 ) 
in Eq.f l44j) . This is because when T — Y 0, the scattering can only happen on the fermi 
surface otherwise it will be Pauli blocked (this is imposed by the vanishing of the prefactor 
f^f^F^Fp). But scattering on the fermi surface yields A a (ki) = A b (k 2 ) = A c (k 3 ) = A d (p) 
and thus K = 0. Therefore, the dimensionless combination K contributes at (D(T 2 /fi 2 ) in 
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FIG. 3: (color online) (rj/s)g 4: In g 1 as functions of (fi/T) 2 for different iVy. 
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FIG. 4: (color online) ((/a 2 T 3 )lng 1 as functions of (fi/T) 2 for different Nf. 



Eq.( l44"j) . which leads to £ oc T. (A similar argument can be applied to rj. As T — >■ 0, f^' = 
fkjL F k 3 F i [ B !j(h) + 4(& 2 ) - B&(fc 3 ) - 5j(p)] 2 = O(T ) in Eq.([M]), since the scattering 
on the fermi surface does not have to be forward scattering thus K' does not have to vanish. 
This leads to r] oc 1/T in Eq.( l36l) .) In contrast to rj, ( is not monotonically increasing with 
Nf because ( is suppressed by an additional power of (T^) 2 oc j3 2 (g 2 ) oc (33 — 2Nf) 2 . Thus 
when Nf is small, ( increases with Nf because quarks tend to make the averaged effective 



13 



S 0.04 



0.02 rr--- - 
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FIG. 5: (color online) (£/ saving 1 as functions of (/jl/T) 2 for different Nf. 

coupling weaker but at large Nf, the suppression factor (3 2 (g 2 ) takes control to make ( 
decrease with Nf. The maximum ( happens when Nf = 5 or 6, depending on the value of 
Ai/T. 

Around fi = 0, we make a Taylor expansion C = a c + ^c(/V^) 2 + ' ' ' ; where is £ at 
zero quark chemical potential. The values of and for varies Nf is tabulated in Table 
IIII1 Our ac is identical to ADM's at least to the second decimal place for all Nf computed 



in Ref. 



17j |. The agreement is better than 1%. 



In Fig. [5j if the coupling g is fixed, our (^/s is monotonically decreasing with Nf. This 
pattern is qualitatively consistent with the massless pion gas result of Ref. [h]] which has 
C/s oc T 4 /N 2 f* oc l/N 2 N 2 (only valid for N f > 2 where pions exist). 



III. CONCLUSION 



We have calculated the shear and bulk viscosities of a weakly coupled quark gluon plasma 
at the leading-log order with finite temperature T and quark chemical potential /i. We have 
found that when normalized by the entropy density s, r]/s increases monotonically with /i 
and eventually scales as (fi/T) 2 at large \i. However (/s is insensitive to \i. Both rj/s and 
(/s are monotonically decreasing function of the quark flavor number Nf when Nf > 2. 
The same property is also observed in pion gas calculations. Our perturbative calculation 
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suggests that QCD becomes the most perfect (with the smallest 77/s) at // = and Nf = 16 
(the maximum Nf with asymptotic freedom). It would be interesting to test whether the 
currently smallest rj/s computed close to the phase transition at /i — and Nf = can be 
further reduced by increasing Nf. 
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TABLE I: Matrix elements squared for two particle scattering processes in QCD. The helicities and 
colors of all initial and final state particles are summed over. q\ and qi represent quarks of distinct 
flavors, qi and q\ are the associated antiquarks, and g represents a gluon. dp and d,A denote the 
dimensions of the fundamental and adjoint representations of SU C (N) gauge group while Cf and 
Ca are the corresponding quadratic Casimirs. In a SU C (3) theory with fundamental representation 
fermions, dp = Ca = 3, Cp = 4/3, and d>A = 8. 



ab — >■ cd 



9i92 -> 9i92 



9i92 ->■ 9l92 



qiq2 ->■ qm 



qiq2 ->■ qiq2 



qxqi ->• 9191 
9i 9i ->• 9i9i 



9191 -> qiqi 



qiqi -> 9292 



9i qi gg 



qw -> qig 
qig -> 919 



95 -> 99 



M, 



a(fei)6(fc2)-5-c(fc3)d(fc4) I 



j2 ^2 

d A 



s 2 + u 2 



d 2 C l( s 2+ u 2 s 2 + t 2 



f ^p 

d A 



t 2 



+ 



l^dpCp (C F - C A /2) 



tu 



4 d 2 F C F ( s 2 + u 2 . t 2 + u 2 



d A V t 2 



+ 



+ l6g 4 d F CF (C F - C A /2) 



u 



,d 2 F C 2 F t 2 + u 2 

*~A ~ 

OA §1 



8g 4 d F C 2 F ( f + ~) - Sg^dpCpCA . ., 



t 2 + u 2 



fdpC 2 (- + -)+ Sg^dpCpCA 



s u 



s 2 + u 2 

t 2 



, . su st tu 
l6g*d A C A [3- ¥ -^-^ 



Appendix A: Scattering Amplitudes and Taylor Expansion Coefficients of Viscosi- 
ties 



To describe the microscope scattering processes in a quark gluon plasma, we need scat- 
tering amplitudes between quarks and gluons. In a hot QCD plasma the infrared singularity 
in the amplitude can be regularized by hard thermal loop dressed propagator. We use the 
same amplitude (shown in Table [I]) as Ref. 
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TABLE II: First two coefficients in the Taylor expansion f\ = + b^n/T) 2 + 
result is identical to AMY's fl^ ] to at least the second decimal place. 



near [i = 0. Our 



N f 
iy i 















27.125 





9 


172.564 


50.381 


1 


60 808 


16.619 


10 


178 839 


51.301 


2 


86.472 


27.281 


11 


184 389 


52.028 


3 


106.664 


34.454 


12 


189.333 


52.608 


4 


122.957 


39.459 


13 


193.764 


53.074 


5 


136.380 


43.055 


14 


197.760 


53.450 


6 


147.627 


45.703 


15 


201.380 


53.755 


7 


157.187 


47.690 


16 


204.675 


54.003 


8 


165.412 


49.207 









TABLE III: power expansion coefficients of Q where Q = + 6^(/i/T) 2 + • • • . Our result is identical 



to ADM's result 



13] to at least the second decimal place. 



N f 




k 


Nf 


a c 


k 





0.4430 





9 


0.3847 


0.0873 


1 


0.5816 


0.0393 


10 


0.3113 


0.0725 


2 


0.6379 


0.0747 


11 


0.2389 


0.0569 


3 


0.6568 


0.0981 


12 


0.1706 


0.0414 


4 


0.6495 


0.1116 


13 


0.1096 


0.0270 


5 


0.6218 


0.1172 


14 


0.0592 


0.0148 


6 


0.5778 


0.1163 


15 


0.0225 


0.0057 


7 


0.5213 


0.1103 


16 


0.0026 


0.0007 


8 


0.4558 


0.1003 









Around \i = 0, we make Taylor expansions fj = (r//T 3 )g 4 lng -1 = a v + b v (/i/T) 2 + - • • 
and ( = (C/a 2 T 3 ) lng -1 = + (yu/T) 2 + • • • . The values of a v , b v , and b^ for various 
Nf are tabulated in Table [Til and IIII1 Our a v and are identical to AMY's and ADM's to 
at least the second decimal place for all Nf computed in Ref. Q and 0, respectively. The 



17 



agreement is better than 1%. 
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